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We show that the isotropic 3-wave kinetic equation is equivalent to the mean field rate equa- 
tions for an aggregation-fragmentation problem with an unusual fragmentation mechanism. This 
analogy is used to write the theory of 3-wave turbulence almost entirely in terms of a single scal- 
ing parameter. A new numerical method for solving the kinetic equation over a large range of 
frequencies is developed by extending Lee's method for solving aggregation equations. The new 
algorithm is validated against some analytic calculations of the Kolmogorov-Zakliarov constant for 
some families of model interaction coefficients. The algorithm is then applied to study some wave 
turbulence problems in which the finiteness of the dissipation scale is an essential feature. Firstly, 
it is shown that for finite capacity cascades, the dissipation of energy becomes independent of the 
cut-off frequency as this cut-off is taken to infinity. This is an explicit indication of the presence of 
a dissipative anomaly. Secondly, a preliminary numerical study is presented of the so-called bottle- 
neck effect in a wave turbulence context. It is found that the structure of the bottleneck depends 
non-trivially on the interaction coefficient. Finally some results are presented on the complemen- 
tary phenomenon of thermalisation in closed wave systems which demonstrates explicitly for the 
first time the existence of so-called mixed solutions of the kinetic equation which exhibit aspects of 
both Kolmogorov-Zakharov and equilibrium equipartition spectra. 

PACS numbers: 



I. INTRODUCTION 

The 3-wavc kinetic equation is the analogue of the 
Boltzmann equation for an ensemble of nonlinear disper- 
sive waves interacting weakly via a quadratic nonlinearity 
in the wave equation. Such wave systems are often con- 
veniently modeled using Hamiltonian equations for the 
complex wave amplitudes, supplemented with addi- 
tional terms modeling forcing, /k, and dissipation, Ft: 



spectral wave-action density, rik, which, for statistically 
homogeneous wave fields, is obtained from the two-point 
correlation function of the wave amplitudes: 



(akiflks) = nki '^(ki - k2). 
It takes the form 



dt 



S[nk] + F[nk] - D[nk]. 



(4) 
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The Hamiltonian, H, contains quadratic and cubic terms 
in the wave amplitudes, a\^: 



H ■ 



I 



j u(k)rfk. 



(2) 



where 



U{]s.i)=j ikikaka (akiakaCtka + akiOkaOks) <5(ki-k2-k3) dk2 

(3) 

The theory of weak wave turbulence [1, 2] studies the 
statistical evolution of the solutions of Eq. (1) in the 
situation where the nonlinear term can be treated as a 
perturbation. Within the framework of weak wave tur- 
bulence, the 3-wave kinetic equation is derived as a con- 
sistent asymptotic closure of the cumulant hierarchy gen- 
erated by Eq. (1). It describes the time evolution of the 
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F[nk]j which is absent for decay problems, represents the 
wave forcing. D[nk] represents dissipation and is typi- 
cally only present at high wave- vectors, ^[nk], referred 
to as the collision integral, describes the conservative 
transfer of energy between wave modes due to resonant 
interactions and is the term responsible for the energy 
cascade. We shall study its explicit form in Sec. II. 

It is convenient, and physically relevant, to consider 
isotropic scale invariant systems. This means that both 
the dispersion relation and the nonlinear interaction co- 
efficient are homogeneous functions of their arguments 
cSigthout any preferred direction. We denote their degrees 
of homogeneity by a and (3 respectively: 



Wafc 
-^akiak2ak3 
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For convenience, we shall take the dispersion relation to 
be a simple power law: 



(7) 



A huge amount is known about the stationary solu- 
tions of Eq. (5) in the turbulent regime where the forc- 
ing and dissipation scales are asymptotically separated 
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from each other in scale. In addition to the thermody- 
namic equilibrium solution, there exists an exact station- 
ary non-equilibrium solution, known as the Kolmogorov- 
Zakharov spectrum, which carries a constant flux of en- 
ergy through scales. This energy cascade solution is the 
analogue of the (phenomenological) Kolmogorov k~^^^ 
spectrum of hydrodynamic turbulence. The fact that 
the energy cascade spectrum can be derived analytically 
is one of the principle reasons for theoretical interest in 
weak wave turbulence. 

Rather less is known about the solutions of Eq. (5) 
beyond the characterisation of the stationary state in 
the limit where the forcing and dissipation scales tend 
to zero and infinity respectively. In particular, knowl- 
edge about the dynamical evolution of the solutions is 
restricted to a subset of systems for which a self-similar 
solution can be constructed using energy conservation 
arguments which in any case, leave the scaling function 
undetermined. Similarly relatively little is known about 
how the system matches itself to the source and sink 
in the case of finite forcing and dissipation scales which 
break the scale invariance necessary to obtain the K-Z 
solution. It is in this context that the present work fits. 

This article contains two main ideas. The first is 
that in the case of isotropic systems, the 3-wave ki- 
netic equation is equivalent to the rate equations for 
an aggregation-fragmentation problem with a rather un- 
usual fragmentation process. This is useful for several 
reasons. Firstly, there is a large body of knowledge about 
rate equations for aggregation-fragmentation problems 
which might provide useful insights. Secondly, this de- 
scription is very compact with almost all properties of 
the solution being determined by a single scaling param- 
eter. Thirdly, and this forms the basis for the second 
main idea of the article, this description forms the ba- 
sis for a new numerical procedure for solving the 3-wave 
kinetic equation which can resolve a very large range of 
scales compared to a direct numerical integration. This 
numerical procedure can then be used to investigate as- 
pects of 3-wave turbulence which are less amenable to 
analytic understanding. 

The layout of the article is as follows. We first ex- 
plain in Sec. lithe analogy between 3-wave turbulence 
and aggregation-fragmentation equations, with the tech- 
nical details relegated to an appendix. One of the prin- 
ciple insights provided by this analogy is that most of 
the scaling properties of the system are determined by a 
single parameter, A = . In Sec. Ill we devote some 
time to explaining how the standard results of weak wave 
turbulence are expressed in terms of A. In Sec. IV we con- 
sider truncating the system at some finite frequency and 
discuss open and closed truncations, two distinct natu- 
ral choices of truncation. In Sec. V the results of Sec. II 
and Sec. IV are used to develop a new numerical pro- 
cedure which allows the stable integration of the 3-wave 
kinetic equation over many decades of frequencies. Some 
technical details of the method are postponed to a sec- 
ond appendix. The remainder of the article is then de- 



voted to presenting some preliminary studies which are 
intended to demonstrate the usefulness of this algorithm. 
In Sec. VII we study the numerical signature of the dis- 
sipative anomaly in finite capacity cascades. In Sec. VIII 
we demonstrate the non-trivial structure of the bottle- 
neck effect in the 3-wave kinetic equation with an open 
truncation and in Sec. IX we study the thermalisation 
phenomenon which occurs when the equation is subjected 
to a closed truncation. The article closes with some con- 
clusions and speculations about future directions of re- 
search. 



II. FORMULATION OF THE ISOTROPIC 
3-WAVE KINETIC EQUATION AS AN 
AGGREGATION-FRAGMENTATION PROBLEM 

The collision integral is usually written in the form 

S[n\^ = I (iJkikskg - -Rkjkaki - -Rkgkika) d^2dk2 



(8) 



where 



^kikaka = 47rLkik2k3('^k2"k3 - "-kinka - nkinks) 

^(wki - Wk2 - Wks) '^(ki - k2 - ks) 

The total wave-action, N, and total quadratic energy, E, 
in the system are 



N 
E 



= I riki dki 
= / Wki nki 



(9) 
(10) 



respectively. E is conserved by Eq. (5) in the absence of 
forcing and dissipation. For isotropic systems it is con- 
venient to work with the angle-averaged frequency spec- 
trum, N^^ instead of the basic k- space spectrum, n\^. 
is defined such that J^^ dui is the total wave action in 
the frequency band [ui , 0)2] ■ It is shown in the appendix 
that, for isotropic systems, Eq. (5) is equivalent to 



dt 

where 



= Si[N^]+S2[N^]+S3[N^]+F[N^]-D[N^] (11) 



Si [N^^]=jKi{ui2,Ul3) N^^N^J{uJi-Ul2-C03) dc023 

-jKi{u3,u)i) N^^N^^5{w2-u)3-i^i) ^^23(12) 

'Ki (Wl, W2) N^^N^J{UJ3-UJI-UJ2) duJ23, 



S2[Nu;] ^-J K2{UJ2,UJ3) N^^N^J{UJ1-UJ2-U)3) duj23 

+jK2{u)3,Wi) N^^N^J{u}2-C03-Ui) dtO'^S) 



+ K2{uJi,uj2) N^^N^^5{ijj3-u>i-u;2) duJ25 
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FIG. 1: Graphical representations of the dynamics encoded 
by Eq. (12), Eq. (13) and Eq. (14). 



and 



23 



'UJi —UJ2) dhJ\ 



23- 



F\NJ\ and D\NJ\ are forcing and dissipation terms whose 
exact forms depend on the problem under consideration. 
In these formulae, as shown in the appendix, Ki{uji,ujj), 
K2{uJi,iOj) and i^s (w^ , ) are homogeneous functions 
which an be constructed from the original interaction co- 
efficient, ikik2k3- They all have degree of homogeneity 



A : 



2/3 -a 



(15) 



We have therefore already obtained a nontrivial result 
which seems to have gone unnoticed before: all scaling 
properties of Eq. (5), which contains 3 scaling parame- 
ters, a,/3 and d, seem to depend on a single scaling param- 
eter, A given by Eq. (15). In fact, this is not completely 
true. Owing to the fact that Ki^ujijUjj), K2{uJi,uJj) and 
K^^LUijLUj) are not identical functions (even though they 
have the same degree of homogeneity) , some memory of 
a,/9 and d is retained in their internal structure - see 
Eq. (105) - which can sometimes be important. The de- 
tails of what can and cannot be expressed solely in terms 
of the new parameter. A, will be addressed in Sec. III. 

Let us first investigate the physical meaning of Eq. (12) 
- Eq. (14). Firstly, it is easily shown that each of the 



three coUision integrals, Si[N^], S2[Nlj] and 5*3 [A''^^], in- 
dividually conserve the total energy. 

Let us first consider Eq. (12) in isolation. The resulting 
kinetic equation 



dt 



(16) 



is actually the Smoluchowski kinetic equation [3] which 
describes the mean-field dynamics of cluster-cluster ag- 
gregation although it has been written in a somewhat 
non-standard form [4, 5]. We can therefore make an anal- 
ogy between wave turbulence and cluster-cluster aggre- 
gation. In this analogy, wave frequency, w, is analogous 
to cluster mass, m, spectral wave-action density, N^^, is 
analogous to density of clusters having mass m, Nm and 
the wave interaction coefficient, Ki{uji^liJ2) is analogous 
to the coagulation kernel, K{mi,m2). 

The physical meaning of the collision integral is easy to 
understand in the aggregation context. In coagulation, 
pairs of clusters having masses nij and aggregate to 
produce a single cluster having mass mrii — nij + mi^. 
The rate at which this process occurs, as for any chemi- 
cal rate equation, is proportional to the aggregation ker- 
nel, K^mj^nik), and to the densities of clusters having 
masses, rrij and rUk- This can be thought of an ex- 
change of mass within the mass triad, (m^, nij, nik). The 
mass contained in clusters of mass rrij decreases at a rate 
mjK{mj,mk)NmjNmki the mass contained in clusters 
of mass mj, decreases at a rate •m]^K{mj,mk)NrnjNmk 
while the mass contained in clusters of mass increases 
at a rate {nij + ruk) K{mj,nik) Nm^. Clearly mass 
is conserved overall. The collision integral therefore sim- 
ply calculates the net rate of change of the density of 
clusters of a given mass by summing the contributions of 
this process over all mass triads. The signs of the terms 
composing the collision integral now make perfect sense: 
the first (positive) term accounts for the rate of increase 
of clusters of mass mi due to the aggregation of pairs of 
smaller clusters having masses 7712 and 7713 which satisfy 
777,2 + "^3 = ™i ■ The second and third (negative) terms 
account for the rate of decrease of clusters of mass nrii 
which occurs when such clusters meet any other cluster 
having mass 7712 or 7*7,3 and aggregate to produce a heav- 
ier cluster having mass ttii -|- m2 of mi + m^. Of course, 
these two negative terms are usually combined into a sin- 
gle term since they differ only by the labeling of dummy 
variables. Nevertheless, for reasons which will become 
apparent in Sec. Ill we choose to keep the rather verbose 
form of Eq. (12). 

Taking this over to the wave analogy, we see that 
the first term on the RHS of Eq. (5) describes a trans- 
fer of energy in resonant triads, (wi, Wj, w^). That is 
to say, triads for which — ujj + uj^.. The energy 
contained in waves of frequency tUj decreases at a rate 
ujjKi{ujj,uJk)N^.N^^, the energy contained in waves of 
frequency uik decreases at a rate uJkKi{ujj,ijJk)Ni^-N^^. 
while the energy contained in waves of frequency uji in- 
creases at a rate (ujj~\-LJk) Ki^LUj^Uk) N^. N^^^- The rates 
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of energy transfer for each mode in the triad {oJi,Wj,oJk) 
can thus be summarised as 

oJk ■ = -uJkKi{wj,Wk)N^jN^^ (17) 

Ui : AEi = +iAEj + AEk). 

This process is illustrated graphically in Fig. 1(A). The 
collision integral then c;alculates the net rate of change 
of the energy of the waves of each frequency by sum- 
ming the contributions of this process over all resonant 
triads. It is clear from this discussion that energy can 
only be transferred from lower frequencies to higher fre- 
quencies by the collision integral Si[Ni^]. As we shall 
see below, the reverse is the case for the integrals 5*2 [A'^] 
and 5*3 [A^^]. <S'i[A^aj] is therefore the driver of the direct 
cascade in wave turbulence. 

If we have laboured the point a little on the inter- 
pretation of 5i[A'ij] and the analogy with aggregation, 
it is because the corresponding interpretation of S2[N^j] 
and iSslAj^] is less obvious. Indeed, the dynamics en- 
coded by these integrals illustrates clearly why wave res- 
onances are different from particles. We first notice that 
the sign structure of 5'2[A^a,] and S^lNi^] is difi'erent from 
that of S'i[A''„]. Two terms are positive and one is neg- 
ative. Looking at the negative term, it is clear that the 
frequency, oji, which is losing energy is the sum of two 
lower frequencies, ui2 and ^3. Thus the correct pictures 
to draw for these processes are those in Fig. 1(B) and 
Fig. 1(C): a high frequency mode loses energy to a pair of 
lower frequency ones, the reverse of the process encoded 
by S'i[At^]. The difference between the two processes is 
in the rates and is summarised as follows. For each reso- 
nant triad, (uji, ujj, Uk) satisfying Wj = uij -|- Wfe, in <S'2[A'j^] 
the rates of energy transfer are: 

Uj : AEj = +Ljj K2 {wj , Wk ) A'cjj N^^^ 

Uk : AEk = +oJkK2{wj,Wk)N^^N^, (18) 

LOi : AE, = -{AEj + AEk) 

On the other hand, for each resonant triad, (wj, w^, w;;) 
satisfying = coj + Uk, in iS'3[A'(^] the rates of energy 
transfer are: 

^3 ■ = +ujjK3{wj,Wk)N^^N^, 
Uk ■■ AEk = +uJkK3{wj,Wk)N^^N^, (19) 
: AE, = -{AE, + AEk) 

It is clear from Eqs. (18) and Eqs. (19) that energy can 
only be transferred from higher frequencies to lower fre- 
quencies by the colHsion integrals S'2[A^^] and S'3[A^:j]. 
Thus they describe back-scatter in the wave turbulent 
cascade. In the aggregation analogy, they can be thought 
of as describing some kind of nonlinear fragmentation 
process in which the rate of fragmentation of clusters of 
a given size is proportional to the density of clusters of 
that size and to the density of fragments. This latter de- 
pendence makes this a rather unusual process from the 



point of view of interacting particle systems. Fragmen- 
tation is often modeled as a linear process. For a review 
of fragmentation see [6] . Although some non- linear mod- 
els have been studied (see, for example, [7] and the ref- 
erences therein) the fragmentation mechanism at work 
here is, to the best of our knowledge, new. The idea that 
nonlinear fragmentation has some connection to energy 
transfer in turbulence is not a new idea [8, 9] but this is 
the first case for which the fragmentation equations can 
be derived from the underlying dynamical equation. 

In what follows, we shall typically work with simpli- 
fied model interaction kernel rather than the complicated 
functions which would arise from particular examples of 
physical wave systems. We introduce the following model 
kernel, which has been very extensively studied [10, 11] 
in the context of cluster-cluster aggregation: 

Ki{uJuC02) = \ Hu^^2 + ^1^^) ■ (20) 

Here the exponents and v must satisfy fi + u = \. Two 
special cases are of particular interest. The first is the 
product kernel: 

Ki{ljJi,UJ2) = {ljJlOJ2)^ ■ (21) 

The second is the sum kernel: 

K^{L0i,i02) = \{uji+uj^). (22) 

In developing the analogy between 3-wave turbulence 
and aggregation-fragmentation problems, it is worth 
pointing out that, in the context if aggregation, it is 
common to work with a discrete analogue of Eq. (16). 
This is for relevant for the case of so-called mono-disperse 
initial conditions which means that all particles initially 
have the same mass which can be taken equal to one. 
The dynamics will only produce clusters having integer 
masses and Eq. (16) can then be presented as an infi- 
nite set of coupled ordinary differential equations for the 
mass densities, Aj, with the integrals having been re- 
placed by sums. A discrete point of view is relevant for 
the wave kinetic equation too. In decay problems with 
monochromatic initial conditions or forced problems with 
monochromatic forcing, it is clear from Eq. (11) that only 
multiples of the input frequency, wq, can be excited by 
the dynamics. Taking = 1? the discrete frequencies 
uji, are then integers, and replacing the integrals with 
sums, we obtain the following infinite set of coupled ordi- 
nary differential equations describing the time evolution 
of the infinite vector of discrete wave occupation num- 
bers, N(t) = (iVi(t),iV2(i),...): 

dN- 

^ = 5i[N]+J-[N]-A[N] i = l,2,... (23) 
where 
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with the discrete colUsion integrals given by 



- Wj-i (24) 



- Y,K^{i,3)NiNj 



i-l 



+ Mj-i,i)NjNj_i (25) 



j=i+i 

oo 



and 



[N]= - ^K3(j,»-i)iV,iV, 



^ if3(j-«,»)^,iVi (26) 



j=i+i 

oo 

+ 5]if3(»,j)^,^^+, 

The forcing and dissipation terms should be chosen ap- 
propriately according to the problem under study. In 
this section we have shown that it is possible to think of 
3-wave turbulence from the perspective of aggregation- 
fragmentation problems. In the remainder of the article, 
we shall demonstrate the usefulness of this analogy for 
understanding turbulence. The analogy turns out to be 
useful in both directions however. For a discussion of 
stochastic aggregation from the point of view of turbu- 
lence theory, see [12, 13]. 



III. A REVIEW OF SOME STANDARD WAVE 
TURBULENCE RESULTS 

In this section we write down some of the standard 
results of wave turbulence theory in the language of the 
previous section. 



A. Kolmogorov-Zakharov Spectrum 

One of the principal results of wave turbulence theory 

is the fact that, for scale invariant systems, the kinetic 
equation has an exact stationary solution which carries 



a constant flux of energy through scales. This solution, 
known as the Kolmogorov-Zakharov (KZ) spectrum, is 
the direct analogue for waves of the well known Kol- 
mogorov k~^^^ spectrum characteristic of hydrodynamic 
turbulence. For isotropic systems, the KZ spectrum is 
usually presented as a stationary solution of Eq. (5) in 
the limit where the forcing wave-number tends to zero 
and the dissipation wave-number tends to infinity. It is 
usually written [1, 2]: 



(27) 



where ^ dimensional constant which can be calcu- 

lated. 

The key step in obtaining the KZ spectrum from 
Eq. (5) is to apply a change of variables known as the Za- 
kharov transformation to the second and third integrals 
in Eq. (5) . The idea of the transformation is to map the 
supports of the frequency delta functions in the second 
and third integrals onto that of the first which allows 
the stationary solution to be clearly seen. The Zakharov 
transformation can be applied individually to each of the 
collision integrals in Eq. (11). We shall demonstrate the 
procedure explicitly for [N^^] and then write down the 
analogous results for S2[N^] and S'3[A''tj]. 

Let us seek a solution of the form iV^^ = ckz with 
X to be determined so that 

-Ckz jKi{ujs,uji) {wiuj^y 5{uj2-oJs-i^i) dhJ2d. (28) 
-Ckz jKi{ijJi,uj2) {ijJiuj2)~'' 5{ijJi-uJi-W2) dhJ2s- 



We now apply the following changes of variables: 

^012' UJ2 ) 



(^^2,^3) 



and 



(^^2,^3) 



( 



UJlljJ2 U) 



W3 



(29) 



(30) 



to the second and third integrals in Eq. (28) respec- 
tively. Noting that the respective Jacobians are 

and and utilising the fact that Ki{uj2,uj3) is a ho- 

mogeneous function of degree A, some algebra yields the 
following: 



0=CkZ JKi{u12,U!3){u}2I^3) ^^(wi - ^2 - ^s) C^W23 



-Ckz Jki{ui2,uj3) ^ 



A+2-2X 



(W2W3) " 

5{ui—U2—C03) duJ23 

X+2-2X 



(31) 



-Ckz / ^1(^2,^3) — (W2W3) ^ 



W3 J 



S{u)i—U)2—C03) dc023- 



6 



This can be put into a single integral, 
0=Ckz jKi{ijj2,uJi) (w2W3)~^^(a;i-a;2-W3) 



\+2-2x r, ,2x-A-2 



— CO. 



2X-X-2 



— Wo 



.2X-X-2- 



(32) 



from which it is easy to see that the right hand side van- 
ishes when 2x — A — 2 = 1. This yields the KZ exponent 
X = The stationary angle-averaged frequency spec- 
trum is therefore 



Noj = ckz w 2 



(33) 



Using Eq. (89) and Eq. (7), it is clear that Eq. (33) is 
equivalent to the more usual expression for the KZ spec- 
trum given in Eq. (27). As is often remarked, Eq. (33) 
can also be obtained simply by dimensional analysis. The 
true worth of the Zakharov transformations lies in the 
fact that they provide a means to obtain the numeri- 
cal value of the constant Ckz and to study the conditions 
under which the spectrum given by Eq. (33) is an admiss- 
able stationary solution of the kinetic equation. Finally, 
returning to the analogy with aggregation, Eq. (33) is 
also well known [4, 14-16] in the aggregation literature 
as the stationary solution of the Smoluchowski equation 
in the presence of a source of monomers. 

Exactly the same steps may be applied to 5'2[A^a;] and 
S'sfA'';^]. The resulting integrals, analogous to Eq. (32) 
are 







K2{(J3,i^2){0JlU;3) ''6{(Ji-(J2-0J3) (34) 



,X+2-2x r, ,2x-X-2 

1^1 



.2X-X-2 



,22:-A-21 



dc023, 



and 



0=-Ckz /■^3(w3,W2) (wia;2) ^(5(wi-u;2-W3) (35) 



,X+2-2x r, ,2x-X-2 



.2X-X-2 



.2X-X-2 



] du)23 



respectively. 



B. Finite and Infinite Capacity Cascades 

The direct energy cascade has infinite capacity [1, 17] if 
the energy contained in the KZ spectrum diverges at the 
high-fc end and finite capacity otherwise. The notion of 
capacity is important because for finite capacity systems 
forced with a constant energy injection rate, the cascade 
necessarily propagates to A; = oo in finite time [18]. In 
the usual notation, the KZ spectrum Eq. (27) has finite 
capacity when j3 > a. In the notation of Sec. II, the 
capacity criterion is be determined by considering the 
energy contained in the KZ-spectrum, Eq. (33) in the 
range of frequencies [wq, f2]: 



/ 



O i-Q, _ A+i 

LoN^^du) = Ckz I ^ ^ w 

0)0 JoJo 



2CKZ 

1- A 



■ Wo 



(36) 
(37) 



Looking at what happens as 17 — > oo. we see that the 
cascade has finite capacity for A > 1. In the aggregation 
analogy, the criterion A > 1 is known as the condition for 
the presence of a gelation transition in the system. 



C. Breakdown Criterion and the Generalised 
Phillips Spectrum 

The derivation of Eq. (5) requires that the linear 
timescale, tl, associated with the waves is much faster 

than the nonlinear timescale, tnl, associated with reso- 
nant energy transfer between waves. This condition may 
be invalidated by the KZ spectrum, either at large or 
small scales [1, 17, 19], a situation referred to as "break- 
down" . The breakdown criterion is derived as follows. 
The linear timescale can be estimated as tl ~ uj^^. The 
nonlinear timescale can be estimated as t-^-^ ^ "^^^ . 
Therefore, on an arbitrary spectrum, ~ oj~^, the ra- 
tio tl/ti^l can be estimated from Eq. (11): 



Tnl 



(38) 



If X is the KZ exponent given by Eq. (33), then this ratio 

becomes 

— ~ (39) 
Tnl 



from which we conclude that the KZ spectrum breaks 
down at high frequencies if A > 3. Using Eq. (15) to 
translate this back into the usual notation, we recover 
the usual criterion [1, 17] for breakdown at small scales, 
/3 > 2a. 

The Generalised Phillips Spectrum [20, 21] is the spec- 
trum for which the ratio tl/tnl is independent of the 
scale. It is important since it is a likely candidate to 
replace the KZ spectrum after breakdown occurs [22] . 
From Eq. (38), it is clear that the Generalised Phillips 
Spectrum in the notation of Sec. II is simply 



(40) 



If this spectrum is translated back into the usual notation 
using Eq. (89) and Eq. (7) we obtain rife ~ fc-(2/3-2a-Kd)^ 
which is the analogue for 3- wave interactions of the better 
known formula for the 4-wave case, Hk ~ fc~('>'~"+'') [20, 
21] (the general formula for N-wave interactions is Uk ~ 

^-(27„-2a-|-(JV-2)d)/(JV-2)^_ 



D. Thermodynamic Spectrum 

In the above application of the Zakhaxov transfor- 
mations to the collision integrals, 5i[A't^], 52[A't^] and 
5'3[Afa;]) we picked out the KZ spectrum as a stationary 
solution in each case. We saw no sign of the other sta- 
tionary solution of Eq. (5) , the thermodynamic spectrum 
corresponding to equipartition of energy. This is because 
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the equilibrium spectrum satisfies detailed balance in the Si[Ni^], S2[N^] or 6*3 [A'j^.] individually but rather of the 
sense that the forward and backward transfer terms bal- full collision integral. To see this, let us add together 
ance each other scale by scale. It is not a property of Eq. (32), Eq. (34) and Eq. (35). The result is 



= CkZ / [^l('^3,t^2) (^2^3) - K2{uJ3,i^2) i^l^^s) - i^3(^'3, ^2) (^1^2)" 



A+2-2x r, ,2x-A-2 



— u. 



,2x-X-2 



2a;-A-21 



(5 (wi — 0^2 — W3 ) (ia;23 



(a;2W3) — ] {^\^^) 

i02 J 



(W1W2) 



,A+2-2x r, ,2x-A-2 



,2x-A-2 



,A+2-2x r, ,2x-A-2 



.2X-A-2 



,2x-A-21 



,2x-A-21 



6{ijJi—uj2—(jJ3) doj- 
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(41) 



(42) 



(43) 



(5(a;i— 0^2— W3) dco23- 



In these manipulations we have used Eqs. (105) and the 
fact that in the integrand. lJ2 + ^^3 = i^'i- It is now 
clear that the total collision integral also vanishes when 
a; = + 1 so that 



(44) 



is also a stationary solution of Eq. (11). Using Eq. (89), 
this spectrum translates into rifc ~ The energy per 
mode is then uj-^nk = const. Hence Eq. (44) corresponds 
to the equilibrium solution. We note that the thermo- 
dynamic spectrum is one aspect of Eq. (5) which cannot 
be expressed in terms of the parameter A introduced in 
Sec. II. 



E. Locality of the Kolmogorov-Zakharov Spectrum 

The Zakharov transformation used to obtain the sta- 
tionary Kolmogorov-Zakharov spectrum is only a valid 
procedure if the collision integral is convergent on the KZ 
spectrum. This property should be checked a-posteriori 
and is referred to as "locality" . The choice of terminology 
comes from the requirement that the collision integral in 
the inertial range should not be dominated by the high 
or low frequency cut-offs. The determination of locality 
is quite a delicate issue so in this section we shall perform 
the analysis exphcitly for Eq. (11). 

We can write Eq. (11) in the form 



dt 



where 



-i?(a;2,a;3,u;i) - ^(^3,c.'i,c^2) C J ' (45) 



R{uJl,UJ2,L0z) = Ki{l03,L02)N^^N^^ 

-K2{uj3,W2)N^,N^,-K3{uJ3,^2)N^,N^^. (46) 



Using Eqs. (105), this can be written 

i?(wi,W2,W3) = Ki{uJ3,U}2) [N^^N^^ 



^ UJi LJ3 



.(47) 



Using the delta functions and integrating out u)3, 
the collision integral can then be written S[N^^] = 
J dco2 1 {1^1,1^2) where 

I{lJi,UJ2) = Q{uJ2,i0i — UJ2) — QiuJ2 — IjOi,U!i) 

-Q{iJi,(j2), (48) 



and 



Q{uji,ijjj) = R{ijJi + Wj,ijJi,Wj). 



(49) 



We need to determine the convergence properties of 
Eq. (48) as L02 ^ Q and UJ2 ^ 00 for a general power 
law distribution, ~ w""^ and then show that the inte- 
gral is convergent when x is the K-Z value. This cannot 
be determined from simply counting powers of uj2 since 
there are some hidden cancellations which occur. Indeed, 
if these cancellations did not occur, so that power count- 
ing would work, it would be impossible to obtain a con- 
vergent integral since no power of LO2 can be integrable 
both at and at 00. 

Let us first examine the limit ci;2 0. We need to 
determine the smallest power of u>2 in I{uJi,u)2) as L02 
0. In this case, the second term in Eq (48) vanishes and, 
after performing a Taylor expansion for small aj2, we find 



I{u;i,(j2) 



*o dQ{y.uj2) 
' 0J2 



dy 



-higher order in UI2 (50) 



We obtain an unexpected cancellation which makes the 
smallest power of CO2 larger by 1 than expected from sim- 
ple power counting. Now we need to determine the be- 
haviour of Q{ljJi,u}2) as a;2 ^ when ~ To 
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accomplish this, wc shall need to know the asymptotic 
behaviour of Ki{oji,ujj). Let us introduce exponents jj, 
and which characterise the asymptotics of Ki{tJi,ojj) 
as follows: 



With some care we shall find that an additional cancel- 
lation occurs within the structure of Q{ljJi,u}2)- When 
= uj-^, Eq. (49) and Eq. (47) give: 



Ki{u)i,iJj) uj^LOj for (jji » a;2. 



(51) 



U}^U2 



2x 



(52) 



Two cases therefore arise depending on whether 1 — a; > 
or 1 - a; < From Eq. (44) it is clear that these 

two cases correspond to the exponent, x, being bigger or 
smaller respectively than the thermodynamic exponent. 
In the former case, the small uj2 behaviour of Q{uji,uj2) 
is 072^^^"'^. Putting this together with Eq. (50), we see 
that the small uj2 behaviour of I{ijJi,uj2) is uj2~'^^'^- In 
this case, the condition for convergence of the collision 
integral at small UJ2'^& f — x + 2> —1 which gives 



X <u + 2). 



(53) 



In the latter case, the small W2 behaviour of (5(wi,W2) 
. The small 0^2 behaviour of I{uj\,uj2) is then 



IS (JO2 



A\ic[\ case, the condition for convergence 



of the collision integral at small ijJ2is v + 
which gives 



d—a 



a — d 



+ l<u + 3. 



+ 1 > -1 



(54) 



Let us now examine the limit UJ2 00. We need to 
determine the largest power of U2- For large, LJ2, there is 
no analogous cancellation between the terms in Eq. (48) 
which led us to Eq. (50). In this case, the first term in 
Eq. (48) vanishes and we find 



the large 0^2 behaviour of I(uji,uj2) is lo2~^~^- In this 
case, the condition for convergence of the collision inte- 
gral at large a;2isA* — a; — 1<— 1 which gives simply 

x> fj,. (57) 

In the latter case, the large W2 behaviour of Q{oJi,u)2) 

u!2^ ° The condition for 

convergence of the collision integral at large u)2 is then 
fj, + — 2x < —1 which gives 



and, hence 



If a-d 
x> - [ ^1+ + 1 



(58) 



The conditions for convergence of the collision integral 
are different depending on whether the exponent, x, is 
greater than or less than the thermodynamic exponent 
given by Eq. (44). In the former case, we must satisfy 
Eq. (57) and Eq. (53). Putting these two together, we 
see that a range of x exists which produce a convergent 
collision integral if 



^Ji<v + 2>. 



(59) 



We note that when such a rangx^ exists, the K-Z exponent 
is at the centre of this range and is therefore local. In the 
latter case, we must satisfy Eq. (58) and Eq. (54) which 
together can be rearranged to give 



/(W1,W2) 



-2Q{uji,W2). 



(55) 



There is still, however, a cancellation within Q(lui,lu2). 
We perform an analysis similar to that leading to 
Eq. (52), except we Taylor expand in wi since it is W2 
which is large in this limit. The result is: 



Q{u}i,U>2) 



W2'-"/(wi)+W2° ^('^i) • (56) 



There are again two casesm depending on whether x is 
bigger or smaller than the thermodynamic exponent. In 
the former case, the large 0)2 behaviour of Q{loi,uj2) is 
Wj"^" • Putting this together with Eq. (55), we see that 



x> -{ii+v + 2,) = Xkz- 



(60) 



Hence we conclude that the K-Z spectrum cannot be lo- 
cal if it is shallower than the thermodynamic spectrum. 
Hence the conditions for locality of the K-Z spectrum are 

XKZ > 1+^^ (61) 
a 

fi < 1/ + 3. 

Note that, provided the K-Z spectrum is steeper than 
the thermodynamic one, the condition for locality only 
depends on the properties of Ki{u)i,u)j) and not on a and 
d. 
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F. Physical Examples 

To close this very brief review of wave turbulence, let 
us calculate the value of A from Eq. (15) for some of the 
commonest physical examples of 3-wave turbulence (see 
[23] for a summary) . Capillary waves on deep water, the 
archetypal example of 3-wave turbulence, have q: = |, 
d = 2 and /? = | giving A = 2. Acoustic turbulence has 
a = 1, d = 3 and /? = | giving again A = 2. Quasi-2D 
Alfven wave turbulence has a = 1, d = 2 and /? = 2 
which yields A = 3. 



IV. CHOICES OF SPECTRAL TRUNCATION 

It is of interest to study Eq. (5) in the presence of a fre- 
quency cut-off, which we shall denote by fl. In the trun- 
cated system, all having ui > Cl are taken to be zero. 
In the discrete case, this corresponds to studying a trun- 
cated version of the infinite set of ODEs constituting the 
kinetic equation. This interest may be forced upon us. In 
a numerical setting, for example, we must necessarily dis- 
cretise and truncate in order to develop a computational 
scheme. Numerical practicalities aside, spectral trunca- 
tions of turbulent systems have recently been of consid- 
erable theoretical interest in their own right because of 
the connection between spectral trunction and the phe- 
nomenon of thermalisation. Thermalisation refers to a 
situation in which a turbulent system exhibits a mixture 
of constant flux and equipartition behaviour. We shall 
have more to say about thermalisation in Sec. IX but 
first, let us c;larify some issues related to the implemen- 
tation of the spectral cut-off. 

We shall truncate the discrete kinetic equation, 
Eq. (23). Requiring that Ni = for i > ft does not 
uniquely determine the resulting set of equations. We 
must choose what to do with loss terms in the forward 
transfer integral, Eq. (23), for which i + j > CI. In the 
graphical representation of Fig. 1, the issue is how to 
treat the set of triads for which Nj and Nk are in the 
truncated system but Ni is not. There is no ambiguity 
arising from such triads in the backscatter terms, Eq. (25) 
and Eq. (26). Prom the rates, Eq. (18) and Eq. (19), that 
if LUi > Ct, then Ni = and all resulting rates are zero. 
There is, however, an ambiguity arising from the forward 
transfer term, Eq. (24). From the energy transfer rates, 
Eq. (17), it is clear that having iVi > CI and Ni = does 
not necessarily imply that the rates are zero although 
terms having u>i > CI can only decrease the total number 
of waves. They correspond to transfer of energy across 
the cutoff from the interaction of two waves which are 
themselves, below the cutoff. Such interactions should 
therefore be treated as dissipation terms in the truncated 
system. Let us separate these terms from the rest and 
write the truncated system as follows, ignoring the ex- 
ternal forcing and dissipation terms for now; 



tZiV ' 

^ = 54N,f2]-7Ti[N,f2] i = 1,2,... Cl (62) 
where 

Si [N, Q] = Sl^^ [N,Cl]+ SP [N, Cl] + [N, Cl] (63) 
with the truncated collision integrals given by 

i-l 

Sl'\N,n]= Y.^,{j,i-j)N^Ni_^ 

j=l 

a 

- K,{j-i,i)NiNj_i (64) 
n-i 

- Y,K^{i,j)NiNj, 

i-l 

sf\^,Cl]= - Y,K2{3,i-3)NiN^ 

Q 

+ K2{j-i,i)NjNj_i (65) 

j=i+i 
n-i 

+ Y,K2{h3)NiNi+^ 

and 

i-l 

5f [N,n]= - Y,K^i3,i-j)NiN,_, 

+ Y K3{3-hi)NjNi (66) 
+ ^i^3(*,j)iV,.iVi+,. 

J=l 

The dissipation terms discussed above corresponding to 
transfer of energy across the cutoff have been gathered 
together into 

(n+i 
Y K,{j-i,i)NiNj_i (67) 
j=n+i 

n \ 
+ Yl K,{i,j)NiNj\, 
j=n-i+i J 

We can now select between different natural truncations 
by varying the parameter 7 in Eq. (62). Taking 7 = 

corresponds to discarding triads which transfer energy 
across the cutoff. We shall refer to this as the closed 
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truncation. For the closed truncation, the energy, Eq, of 
the truncated system, 

n 

En = J2i^i' (68) 

i=l 

is conserved by Eq. (62). Taking 7=1 means that we 
allow energy to freely cross the cutoff at which point it is 
removed from the system (dissipated). We shall refer to 
this as the open truncation. With the open truncation, 
the energy, Eq, of the truncated system may decrease 
as a function of time. Taking < 7 < 1 corresponds 
to something intermediate between the open and closed 
truncations which we shall refer to as a partially open 
truncation. 

There is nothing intrinsic to the system which tells us 
which truncation we should choose. It will depend on 
what we want to do. It is often the case that we consider 
the truncated system as being an approximation of the 
original kinetic equation and would like to recover the 
original dynamics when we take ^ 00. As we shall 
see in Sec. VII, the choice of truncation is sometimes 
irrelevant for recovering the original dynamics as — > 00 
and sometimes essential. 



V. A NEW NUMERICAL ALGORITHM FOR 
SOLVING THE 3- WAVE KINETIC EQUATION 

Recasting the 3-wave kinetic equation, Eq. (5) as an 
aggregation- fragmentation problem has another advan- 
tage in addition to the conceptual clarity discussed in 
Sec. III. It is the basis of a new numerical method for 
solving Eq. (5) accurately over a large range of frequency 
scales. This method is an extension to Eq. (11) of an 
elegant method developed by M.H. Lee [24, 25] to solve 
the Smoluchowski equation, Eq. (16). In this section we 
shall give the details of this method. 

Before proceeding, we would like to remark that the 
objective is to design a numerical method which can 
solve the isotropic kinetic equation, Eq. (11), accurately 
over many decades of frequency for an arbitrary interac- 
tion coefficient with the objective of studying the scaling 
properties of this idealised system. This is in contrast 
with, but complementary to, the majority of numerical 
effort invested in solving wave turbulence kinetic equa- 
tions which has focused on approximating the 4- wave ki- 
netic equation for the specific case of deep water grav- 
ity waves under anisotropic conditions [26-28] which is 
strongly motivated by the applications to wave forecast- 
ing. 

Any numerical scheme requires that we work with a 
set of N discrete frequencies, Wj = wq + (« — 1) Aw, where 
i = I, . . . N, Alo is the mode spacing and lon = is the 
frequency cut-off. We shall therefore take loq = I and 
Aw = 1 and work with Eq. (23). It is more natural to 
work with the energies of the modes, rather than their 
corresponding occupation numbers, N. E(i), the vector 



of energies is obtained from N(f) by the relation E^ = 
iNi, i = 1, . . . , n. 

The objective is to solve the following set of coupled 
nonlinear ordinary differential equations for the E^ ob- 
tained by multiplying Eq. (23) by i: 

dE 

= iSi[N,Q]-jiT[N,Q] (69) 
+iFi[N]-iDi[N] 

with ^^[N, fi] and r[N, Q] given by Eq. (63) and Eq. (67) 
respectively. Depending on the application, it might also 
be useful to keep track of the cumulative energy, Eoit), 
dissipated by the T[N, 17] terms. In such cases, we sup- 
plement Eq. (69) with the following equation: 

^=7E*T,[E]. (70) 

We shall take the forcing to be 

Fm = JSi,i (71) 

so that the total rate of injection of energy into the sys- 
tem is J. We take Di[N] to be zero since dissipation is 
provided by T[N, f2]. Then setting 7 = 1 we recover the 
open truncation. Setting 7 = we recover the closed 
truncation and by setting 7 somewhere between, we get 
a partially open truncation. 

There are several problems to be overcome in solv- 
ing Eqs. (69). Firstly Eqs. (69) become very stiff as 
A increases. This means that when straightforward ex- 
plicit integration routines are applied to the system, the 
time step required to maintain numerical stability re- 
mains very small even when the solution we are trying 
to compute is varying slowly. The result is that it takes 
an impractically long time to compute the solution with 
explicit methods and one must resort to an implicit in- 
tegration algorithm in many cases. The second problem 
is that we are interested in situations when the number 
of modes involved is large but the number of modes one 
can deal with by direct integration of Eqs. (69) is limited 
by the computational cost of evaluating the sums on the 
righthand side. This issue is addressed by the technique 
developed in [24]. The modes are grouped into exponen- 
tially spaced bins and the net exchange of energy between 
triads of bins is approximated rather than the exchange 
of energy between triads of individual modes. 

A. Time-stepping 

The solutions of Eq. (69) exhibit scaling behaviour, 
with the result that there is a wide variation of the 

timescale during the course of the evolution. This ne- 
cessitates [24] the use of adaptive timestepping to keep 
the error within a prescribed limit [29]. The explicit 
second order Runge-Kutta (RK2) method and implicit 
trapezoidal (IT) method were both used in conjunction 
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0.01 



Implicit Trapezoidal Method 
■f^ Order Runge-Kutta Method 




FIG. 2: Illustration of the stiffness of Eq. (69) with the A = 2 
sum kernel, Eq. (22). Q, was taken to be 100, the energy injec- 
tion rate was J = 1 and the open truncation was used. The 
plot compares the timesteps required to keep the numerical 
solution within an error tolerance of 1 x 10~® using the RK2 
algorithm and IT method. 



with a step-doubling procedure to adjust the timestep, 
h. Both have stepwise errors of 0{h'^). The use of an 
implicit method is necessary for larger values of A be- 
cause Eq. (69) becomes increasing stiff for larger values 
of A. This is illustrated clearly in Fig. 2 which compares 
the stepsize required by the explicit RK2 and implicit 
IT integration algorithms in order to maintain a given 
error as the solution of Eq. (69) with the sum kernel, 
Eq. (22) with A = 2. Initially, the dynamics is fast an 
both methods require small timesteps to keep the error 
within the prescribed tolerance. Around t — 2 the so- 
lution approaches the stationary state and the dynamics 
slows down. The explicit RK2 algorithm continues to re- 
quire small steps despite the fact that the solution is no 
longer evolving quickly. In contrast the implicit IT algo- 
rithm can take increasingly large steps as the dynamics 
slows down. This behaviour is the classical symptom of 
a stiff system and renders explicit solvers practically use- 
less for solving the 3WKE with larger values of A. 

There is a price to be paid for dealing with the stiff- 
ness issue. The IT method requires that we solve the 
following set of implicit nonlinear equations (we ignore 
the forcing and dissipation terms for now since they are 
straightforward to include) to find the energies at the 
next timestep, Ei(t + h), from the energies at the current 
timestep, Ei(t): 

E,{t + h) = E,{t) + J {S,[^{t)] + 5,[E(t -h h)]) (72) 



This was done using the GSL implementation [30] of the 
Roscnbrock algorithm [31], a standard method of multi- 
dimensional root finding. The current values, Ei{t), were 
used as the initial guess for the root finding procedure. 



B. Coarse-graining 

A direct integration of Eq. (23) is practical only for 
relatively small numbers of modes. In order to resolve 
large inertial ranges, we coarse-grain the modes into bins 
and approximately compute the net energy transfer rate 
between bins following the approach developed in [24] for 
the Smoluchowski equation. 

We need to divide the frequency domain, [l,ri)], up 
into N bins. We shall adopt the notation 

to denote the bin and denote the bin widths by 

We use the same bin structure as adopted in [24]: the 
first n bins are linearly spaced and the next N ^ n bins 
are defined by a geometric sequence of boundary points 
having ratio a — 10". The final bin is defined so that 
its right boundary is fi. With this definition, there are 
approximately n bins per decade of frequency space with 
the total number of bins, N, determined by the value 
of the frequency cut-off, fi. We define a characteristic 
frequency, f^i, of each bin by 

We shall continue to use Ei to denote the total amount of 
energy contained in bin i despite the fact that for i > n 
this now refers to an entire bin of frequencies rather than 
an individual frequency (for the first n bins, each of which 
contains only a single mode, Ei remains just the energy 
contained in that mode). Likewise we shall continue to 
use Ni to denote the total number of waves in bin i. 
Within each bin, Bi, with i > n the energy and wave- 
action distributions within the bin are approximated with 
power law distributions: 



e^(w) = a, ( j < < wf (73) 

n,;(w) = a, ( -^j . (74) 

The exponent, bi, is obtained by interpolating the char- 
acteristic energies of the neighbouring bins: 



log(fe)-log(fe) 



logJ]i+i - logf]j_i 

and the pref actor, a^, is fixed by the normalisation 



duj ei{uj). 



(75) 



(76) 



We have seen in Sec. II that computing the collision 
integral on the RHS of Eq. (23) is equivalent to com- 
puting the appropriate rates of energy exchange given 
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by Eq. (17), Eq. (18) and Eq. (19) for the members of 
each resonant triad and summing over all triads. After 
course-graining, we need to compute the net rates of en- 
ergy transfer resulting from the interactions of all triads 
{cOi,Uj,Uk) for which € and uij € Bj (of course 
we must allow the possibility that Bi=Bj) and sum the 
results over all possible pairs of bins Bi and Bj. The fun- 
damental quantity which we evolve is naturally the total 
amount of energy contained in each bin. 

When considering interactions between a pair of bins, 
Bj and B^, we shall adopt the labeling used in Fig. 1: bin 
j to the left of bin k. In computing the rates of transfer 
of energy generated by these interactions the first step 
is to determine which bin or bins contain the modes in 
resonance with those in Bj and B^. Before describing 
how the subsequent energy redistribution works in detail, 
let us explain the key approximation which will used. 
Consider, for example, the forward-scatter process. The 
total rate of energy transfer to higher frequencies due to 
resonances between modes in Bi and Bf. is given by a 
double integral. 



R 

= / d^j / 



dUk {uJj + Ldk) Kl{ujj,UJk) 

aj{uj)nk{oJk), (77) 



and these rates of gain of energy are distributed (non- 
uniformly) among a set of modes with frequencies lying 
between uj' + ojj^ and cjj^ + uj[^'. The key approximation 
which we make is to treat all modes in Bj (the lower fre- 
quency narrower bin) as having frequency ^Ij (of course, 
this is not even an approximation if j <= n). That is 
we replace nj{LOj) with Nj 5{u>j — Qj) so that instead of 
computing double integrals, we need to compute one di- 
mensional integrals of the form 



dujk (Oj -I- Wfc) Ki{nj,u)k) nk{ijJk) (78) 



with the gains of energy distributed among a set of modes 
with frequencies lying between Qj + ojj: and Qj + u)^ . 

With this approximation in mind, let us now describe 
explicitly the calculation of the rates of energy transfer. 
For each pair of bins, Bj and Bk, there are three possi- 
bilities which must be treated separately. They are, as 
listed below, labeled A, B and C . With apologies for the 
somewhat clumsy notation, we now list the energy trans- 
fer rates to/from each bin for each of these cases split up 
according to the contributions from each of the processes 
described by S\[NJ\, 52 [A^^] and Sz[N^]. L shall denote 
an energy loss term and G an energy gain term. The 
actual integrals are written explicitly in appendix XI to 
avoid an overload of technical details. 



• Case A : k < n 
This is the simplest case in which Bk (and hence 

Bj) is a discrete bin. Given j and k, we identify 
the bin, B^, which contains = flj + il,k- Whether 



it is a discrete or continuous bin is irrelevant. The 
transfer rates can be read off almost immediately 
from Eqs. (17), Eqs. (18) and Eqs. (19): 



(Ai?), 

iAE)k 
(Ai?), 



;(2) 

.(2) 
'k 

(AE)^ + {AE)k 



(3) 

i 

(3) 



(79) 



where 



Aril) 

^] 
Aril) 

AQi2) 



njKi{nj,nk)NjNk 
nkKi{nj,nk)NjNk 



'G 



(2) 



(3) 



^G^ 



Aq{3) 



Q.jK2{Q.j,Q.k)Nj 

nkK2{n,,nk)Nj 
njKs{nj,nk)Nk 
nkK3{nj,nk)Nk 



AtOi 

Ni 
Auji 

A, 
Awi 

Auji 



• Case B : fc > n and \Q.j+w^,Q.j + oj^] is contained 
in a single bin, Bi 

In this case, we identify the bin, Bi, which contains 
all modes in the range ^j+w^, fij-t-Wj^]. The rates 

of energy transfer are 



{AE), 

{AE)k 
{AE), 



{AE)j 



{AE)k 



(3) 

j 

(3) 



(80) 



with the relevant integrals provided in the ap- 
pendix. 



• Case C : k > n and \p.j 



'^k ' 



j^] is split 



between two bins, and Bi^ 
It can happen that the range of modes [fij -I- 
wf', fij -I- u^] spans two bins. With the above def- 
inition of the bin structure it is never more than 



two. We identify these bins as Bi^ and Bi 



We 



must also identify the point Q.* = tof^ — Qj G Bk 
which marks the boundary between those modes 
in Bk which are resonant with modes in Bi^ and 
those modes in Bk which are resonant with modes 
in Bij^ . Energy transfer is then split appropriately 



between Bij^ and 



(AE), 
{AE)k 
(AE),, 
{AE)i, 

where 

L{AE)j = 



l{AE),+r {AE), 
L{AE)k+R {AE)k 
l{AE)j+l {AE)k 
R{AE)j +R {AE)k 



(81) 



Cr(l) 



C ^1,. 



(2) 



C ^1,. 



(3) 
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FIG. 3: Sample output from implementation the algorithm 
described in Sec. V. This is a decay simulation with A = 3/2 
product kernel, Eq. (21). 



with the measurement of stationary scaling exponents, 
the measurement of ckz is a stringent test which the 
code was required to pass. 

Let us now calculate ckz- Consider the statement of 
conservation of energy. 



(83) 



where is the energy flux at frequency lo. For a general 
power law spectrum, = co;"^, we can write the right 
hand side like Eq. (41). Introducing rescaled integration 
variables, u and v defined by 

UJ2 — ^iV OJ3 = UJi U, 

and integrating out v allows us to write Eq. (83) as 



8 J — 



2 ,A-2x+2 



I{x), 



(84) 



where 



LiAE)k 



Cril) 
Cr{l) 



C ^(2) 
'R 

C ^(2) 

C M2) 
'R '-'k 



C ^(3) 

-R 



(3) 
k 



(82) 



C ^(3) 
'R '-'k ■ 



Again, the actual integrals are written explicitly in 
the appendix. 

Some typical outputs resulting from the implementa- 
tion of the algorithm described in this section is shown 
in Fig. 3. 



VI. NUMERICAL VALIDATION: 
CALCULATION OF THE 
KOLMOGOROV-ZAKHAROV CONSTANT 

The evolution of the energy spectrum shown in Fig. 3 
certainly looks plausible. Nevertheless, the approxima- 
tions made in deriving the algorithm described in Sec. V 
are not systematic and no convergence results have been 
proven. Therefore it is essential to validate the code. 
From the discussions of Sec. II, it is obviously inadequate 
to use the measurement of scaling exponents as a means 
of validation since the scaling properties of Eq. (11) and 
Eq. (16) are practically identical even though the physics 
is very different for the two equations. For Lee's original 
implementation of this method for Smoluchowski equa- 
tion, one had the luxury of several exact solutions [3, 24] 
against which the full time-dependent evolution could 
be validated. There are, to the best of our knowledge, 
no known exact solutions of the 3-wave kinetic equation 
whch could play a similar diagnostic role in the present 
context. We suggest instead, to use the measurement of 
the Kolmogorov-Zakharov constant, ckz, as a diagnostic. 
It can be computed exactly as we shall now show. Fur- 
thermore its value is dependent on getting the internal 
structure of the collision integral correct. Taken together 



I{x) 



[Ki{uA-u) (it(l -u))-^ 



(85) 



-K2{u, 1 - u) - i^3(u, 1 - u) (1 - uy] 



Integrating once, we get 



,2x-\-2'\ 



du. 



,A-2a;+3 



A - 2.T + 3 ■ 



(86) 



The energy flux, J^, should be constant, independent of 
LO, and equal to J, the rate of energy injection, when x 
takes the KZ value, (A-t-3) /2. The limit needs to be taken 
using I'Hopital's Rule since /(xkz) = 0. The result is 



3 



1 , dl 

- c — 

2 dx 



(87) 



X=£CkZ 

The K-Z constant, is therefore given by 



Ckz 



/2 J 



dl 

dx 



-1 



The integral ^^r^-^^^i can be calculated numerically, 
and analytically for some interaction coefficients. We 
shall use this result to validate the numerical solution 
procedure. 

Based on the discussion of the previous section, very 
few properties of the kinetic equation actually depend on 
a or c? and those that do are all related to thermody- 
namic equilibrium. We shall therefore take a = c? so that 
Ki{uji,ujj) = K2{uji,ujj) = K^{uji,ujj). For such systems, 
the thermodynamic equilibrium spectrum is uj~^ . We 
tested our code with both the product kernel, Eq. (21), 
and the sum kernel, Eq. (22), for various values of A. 
This was done by allowing the system to reach a sta- 
tionary state with a constant rate of injection and then 
compensating the computed spectrum, iV^ , with the K-Z 
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FIG. 4: Numerical values of the Kolmogorov-Zakharov con- 
stant for the sum kernel, Eq. (22) and the product kernel, 
Eq. (21), for a range of values of A. The theoretical curves 
come from numerical evaluation of Eq. (88) . 



FIG. 5: Time dependence of the dissipated energy for increas- 
ing values of the dissipation scale, Q, for the product kernel, 
Eq. (21) with A = |. 



spectrum and fitting a constant to the result. The results 
are shown in Fig. 4. The numerically fitted values of ckz 
are in excellent agreement with the theoretical predic- 
tions obtained from Eq. (88) for both kernels which gives 
us confidence in the numerical algorithm introduced in 
Sec. V. Note that the value of the ckz diverges at A = — 1 
for both kernels. This comes from the fact that, in the 
numerics we took a = c? so the the thermodynamic spec- 
trum, Eq. (44) has exponent x = 1. Thus at A = — 1 the 
thermodynamic and K-Z exponents coincide. From the 
locality conditions, Eq. (61), derived in Sec. Ill, we know 
that the collision integral diverges when this occurs due 
to the violation of the first condition in Eq. (61). No- 
tice also that the vanishing of ckz for the sum kernel at 
A = 3 corresponds to the breakdown of locality due to 
the violation of the second condition in Eq. (61). Fig. 4 
therefore provides a nice consistency check on our earlier 
theoretical analysis. 
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FIG. 6: Time dependence of the dissipated energy for increas- 
ing values of the dissipation scale, Q,, for the product kernel, 
Eq. (21) with A = |. 



VII. FINITE AND INFINITE CAPACITY 
CASCADES AND DISSIPATIVE ANOMALY 

The concept of a dissipative anomaly is key to un- 
derstanding the statistical properties of many turbulent 
systems. It refers to a situation in which the average rate 
of energy dissipation tends to become independent of the 
dissipation parameter in the dynamical equations in the 
limit where this dissipation parameter is taken to zero. 
This phenomenon is explicitly demonstrable in Burgers' 
equation (see [32] for a review) where the rate of dissi- 
pation of energy in shocks becomes independent of the 
viscosity, as v ^ Q. It is also believed to be very rel- 
evant (see the discussion in [33]) for the high Reynolds 
number limit of hydrodynamic turbulence. 

In wave turbulence it is expected that the dissipative 
anomaly is present for finite capacity wave systems and 



absent for infinite capacity systems. In Sec. Ill we showed 
that A = 1 is the boundary between infinite and finite ca- 
pacity. The difference between the two can be explicitly 
demonstrated using the numerical scheme developed in 
Sec. V. We have seen that an open truncation of the 
3 wave kinetic equation introduces dissipation of energy 
with the dissipation scale given by fi, the cut-off fre- 
quency. The non-dissipative limit in this situation corre- 
sponds to taking this dissipative cut-off to infinity. Figs. 
5 and 6 show the total dissipated energy as a function 
of time for two decay systems with the same initial con- 
dition, Nt^(Q) = S{lu — 1), interacting via the product 
kernel, Eq. (21) with two different values of A. Fig. 5 
shows total dissipated energy as a function of time for 
the case A = 3/4 which is infinite capacity. As the cut- 
off, J7, is increased the time at which the energy is dis- 
sipated increases. Consequentially, the dissipated energy 
at any fixed time vanishes as the dissipative cut-off is re- 
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FIG. 7: The bottleneck effect in the 3-wave kinetic equation 
with constant energy injection at a; = 1 and open truncation 
at 57 = 10*'. The plot shows stationary spectra compensated 
by the K-Z scaling for the product kernel, Eq. (21), for a range 
of values of A. The bottleneck effect results in an accumula- 
tion or depletion of the spectrum near the cut-off, depending 
on the choice of kernel. Solid lines indicate the theoretically 
predicted values of the K-Z constant from Eq. 



moved to infinity. Hence there is no dissipative anomaly. 
The corresponding situation for the finite capacity case, 
A — shown in Fig. 6 is strikingly different. As the 
dissipative cut-off, fi, is removed, the dissipated energy 
as a function of time becomes independent of il. From 
the figure it is clear that for times larger than about 3, 
the dissipated energy is finite in the limit w — > oo. This 
is as explicit a demonstration as one may expect from 
a numerical simulation of the presence of a dissipative 
anomaly in this system. 

It is worth noting that the corresponding dissipative 
anomaly for the Smoluchowski kinetic equation is well 
known in the aggregation literature where it is more com- 
monly referred to as the gelation transition. The crite- 
rion, A > 1, is common to both systems. A numerical 
study similar to what has been presented here is provided 
for that system in [5]. An exact solution for the Smolu- 
chowski equation with the product kernel which exhibits 
the anomaly explicitly can be found in Sec. 4.5 of [3]. 
In the light of the numerical results presented here, the 
construction of the corresponding mathematically rigor- 
ous solution of the 3-wave kinetic equation would likely 
be a fruitful line of research. 



VIII. THE BOTTLENECK EFFECT IN 
SYSTEMS WITH OPEN TRUNCATION 

In the previous section, we were interested in the ba- 
haviour of the solution of the 3-wave kinetic equation 
with open truncation at fixed time as the dissipative cut- 
off, ri, tends to infinity. In this section, we shall consider 
the complementary situation where is fixed and time 



tends to infinity. In this situation, irrespective of whether 
there is a dissipative anomaly or not, we always expect 
the solution to tend to a stationary state corresponding 
to the Kolmogorov-Zakharov spectrum, Eq. (33), mod- 
ified by the presence of the finite cut-off. A question of 
interest is how the K-Z spectrum matches to the dissipa- 
tion range and the first issue which arises is whether or 
not there is a "bottleneck" effect. 

The term "bottleneck" refers to a phenomenon, ini- 
tially observed in numerical simulations of Navier-Stokes 
turbulence (see [34] and the references therein) where the 
stationary spectrum exhibits a "bump" super-imposed 
upon the expected constant flux spectrum as it enters 
the dissipation range. A pysical mechanism for the bot- 
tleneck was suggested in [35] . A corresponding theory for 
the wave turbulence case was suggested in [36]. It is an 
intrinsically dissipative phenomenon. It arises because 
the rate of forward transfer of energy from frequency ui 
due to the interaction with frequency W2 > is propor- 
tional to the product, N^-^N^^ , of occupation numbers of 
both frequencies (see Eq. (17)). The uj2 contributing to 
the total rate of transfer of energy can be divided into 
those UJ2 less than the dissipative scale, il, and those 
greater than fl. The latter set are in the dissipative 
range, where N^^^ is effectively zero (exactly zero in the 
case of dissipation via an open truncation as is considered 
in this article) . Therefore the rate of forward transfer of 
energy through lui would be decreased. However since 
uji is still in the inertial range, the rate of energy trans- 
fer must remain equal to the injected flux. Hence the 
occupation numbers of those UJ2 between uji and must 
increase so that the same flux can be carried by fewer 
triads (hence the term "bottleneck"). This effect then 
produces the bump at the end of the spectrum. 

Here we investigate the bottleneck effect in the 3-wave 
kinetic equation explicitly using the numerical algorithm 
developed in Sec. V. Fig. 7 shows some results. We 
computed the stationary state of Eq. (11) with an open 
truncation a,t H. — 10^ and the product kernel, Eq. (21) 
for several different values of A. Fig. 7 shows several such 
stationary spectra compensated by the corresponding K- 
Z spectra. The solid lines indicate the fitted values of the 
K-Z constant used for validation of the code in Sec. VI. 
We see that there is non-trivial structure super-imposed 
upon the K-Z spectrum as it approaches the dissipative 
cut-off. Interestingly, whether this structure corresponds 
to a bump or not depends upon the value of A. The 
heuristic argument outlined above might lead one to ex- 
pect that the bottleneck effect always produces a bump 
whereas, in reality, this does not seem to be the case. 
From the numerical calculations, it seems that the match- 
ing of the spectrum to the dissipation produces a bump 
for A > 1 but produces a hollow for A < 1. We shall 
remain agnostic on the question of whether it is appro- 
priate to call each of these regimes a bottleneck. While 
this is not a systematic study, it suggests that it might 
be worthwhile to revisit the bottleneck phenomenon in 
wave turbulence. 
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IX. THERMALISATION IN SYSTEMS WITH 
CLOSED TRUNCATION 

In the previous two sections, we investigated two phe- 
nomena - the dissipative anomaly and the bottleneck ef- 
fect - which were related to the choice of open truncation. 
In this section we consider the 3-wave kinetic equation 
with closed truncation where a new phenomenon arises 
known as thermalisation. 

We have already seen that Eq. (5) has two station- 
ary scaling solutions - a thermodynamic solution and 
a KZ solution. The former has finite temperature and 
zero energy flux whereas the latter has finite energy flux 
and zero temperature. It is conjectured [2, 37] that the 
general solution of Eq. (5) should be a two-parameter 
function depending on both a flux and a temperature 
although little is known about what such mixed states 
should look like except, perhaps, as pertubations of the 
pure K-Z or pure thermodynamic solutions (see chap. 4 
of [2]). They were first realised in [38] in the context 
of the Leith model - a very simplified model of hydrody- 
namic turbulence - and soon after were properly observed 
in the full Euler equations [39]. Since then, considerable 
effort has been invested in understanding the interplay 
between thermaUsation and turbulence in the hydrody- 
namic context [40 42] where it is argued [40] that the 
closed truncation arises naturally. Despite all of this ac- 
tivity, the phenomenon has not yet been investigated in 
the context of wave turbulence. 

In order for thermalisation to be possible, the energy of 
the truncated system should be conserved. Following the 
discussion of Sec. IV, we should choose 7 = in Eq. (62) 
and work with the closed truncation. We performed a 
series of decay calculations of Eq. (62) with initial con- 
dition N^^ = S{uj — 1) and truncation frequency fl = 10^ 
for the product kernel, Eq. (21), with a range of values 
of A. Some representative results are shown in Fig. 8. 
Since these are decay simulations, the spectra shown in 
Fig. 8 are not stationary. The spectra are presented at 
times for which the corresponding system with an open 
truncation would have dissipated half of the initial en- 
ergy. Thus the results of Fig. 8 are directly comparable 
with the decay simulations used to investigate the dissi- 
pative anomaly in Sec. VII. Fig. 8(A) shows the spectra 
for several different values of A. The solid black lines in- 
dicate the corresponding K-Z spectra. The effect of the 
closed truncation is clearly evident in the increase of the 
spectrum above its K-Z value near the cut-off. That this 
corresponds to thermalisation of the high frequencies is 
evident from Fig. 8(B) which shows the corresponding en- 
ergy spectra. It is clear that that the spectrum is crossing 
over to an equipartition of energy near the cut-off. 

Unlike the dissipative bottleneck associated with the 
open truncation which we studied in Sec. VIII, thermal- 
isation always leads to an accumulation of energy near 
the cut-off, irrespective of the interaction coefficient. The 
distinction is clearly illustrated by comparing the A = 
cases for the closed and open truncations. The closed 



truncation leads to accumulation of energy near the cut- 
off as shown in Fig. 8(A) whereas the open truncation 
leads to depletion of energy near the cut-off as shown in 
Fig. 7. Since thermalisation always leads to accumulation 
of energy near the cut-off, this effect is sometimes also re- 
ferred to as a "bottleneck". In the current study, such 
dual use of terminology would be quite confusing since 
the closed truncation leading to thermalisation has, by 
construction, zero flux at lo = CI. On the other hand, the 
open truncation leading to the bottleneck effect has, by 
construction, a finite fiux at w = fi. 

The two phenomena can be linked to each other using 
the partially open truncation discussed in Sec. IV. Al- 
though it seems unlikely that such a boundary condition 
is of much relevance to any physical system, choosing the 
parameter 7 in Eq. (11) to have intermediate values be- 
tween and 1 allows us to interpolate smoothly between 
the open and closed truncation and hence between small 
scale thermalisation and small scale bottleneck. The re- 
sults of one suc;li exerc;ise in numerical trickery is shown 
in Fig. 9. In this figure, spectra compensated by the 
K-Z scaling are shown which were obtained by solving 
Eq. (11) with constant kernel (A = 0) and several differ- 
ent values of 7 for a system truncated at = 10^. The 
spectra shown for the finite 7 are truly stationary. Un- 
like the spectra in Fig. 8 discussed above, the spectra in 
Fig. 9 were obtained with a forcing term injecting energy 
at a constant rate. Obviously this was not possible for 
the case 7 = where the total energy diverges in the 
forced case. The spectrum for 7 is not stationary and is 
displayed at a time wliic;h allows it to be compared with 
the other spectra. The message to be taken from Fig. 8 
is that the depletion of energy near the cut-off due to the 
open cut-off goes over to a thermalised accumulation of 
energy near the cut-off as the efficiency of the dissipa- 
tion is decreased. It would be interesting to investigate 
the relationship between the open and closed truncations 
more carefully, especially to understand the role played 
by the interaction coeffic-ient in determining the shape of 
the spectrum near the cut-off. 



X. CONCLUSIONS 

To conclude, we have outlined an analogy between 
the isotropic 3-wave kinetic equation and the rate equa- 
tions for a aggregation-fragmentation problem with an 
unusual nonlinear fragmentation mechanism. This anal- 
ogy demonstrates that almost all properties of the system 
are determined by a single scaling parameter, A = 
thus greatly reducing the parameter space of possible 
behaviours. A new numerical scheme was constructed 
based on this analogy which allows for the stable integra- 
tion of the isotropic 3-wave kinetic equation over many 
decades of frequency space. This algorithm was validated 
by comparing numerical measurements of the stationary 
state with theoretical calculations of the K-Z constant for 
a range of model interaction coefficients. Several applica- 
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FIG. 8: Thermalisation in decay simulations of the 3- wave kinetic equation, Eq. (11), with closed truncation [j = 0) at fl — 10® 
for the product kernel, Eq. (21), for several different values of A. The left panel shows the wave spectrum, TV^,, and the right 
panel shows the corresponding energy distribution, E^j, from which the equipartition effect near the cut-off is clear. 




FIG. 9: Partial thermalisation in the 3-wave kinetic equation, 
Eq. (11), with constant kernel and a partially open trunca- 
tion at Q = 10®. The plot shows stationary spectra (with 
the exception of the case 7 = for which there is no station- 
ary state) compensated by the K-Z scaling uj^^^ for several 
different values of 7. 



likely be fruitful. In addition, since the initial studies of 
Galtier et al. [43] of the solutions of the Alfven wave 
kinetic equation, there is growing evidence [25, 38, 44] 
that finite capacity cascade often, if not always, exhibit 
a dynamical scaling anomaly during the transient stage 
of evolution which cannot be understood from elemen- 
tary scaling arguments. The methods developed in this 
article should provide an ideal set of tools to study this 
phenomenon in the setting of general kinetic equations. 
This will form the basis of future work. The other obvi- 
ous line or research which has opened up is the question 
of whether the methods decribed here can be extended 
to the case of the isotropic 4-wave kinetic equation. At 
this point, the answer does not seem obvious. The key 
approximation used in Sec. V to develop the numerical 
algorithm involved treating all waves in the leftmost bin 
as having the same frequency. This may not be a reason- 
able approximation in a 4-wave system where the waves 
interact in quartets and the possibility of an inverse cas- 
cade may increase the sensitivity of the dynamics to the 
way in which low frequencies are approximated. 



tions of the new algorithm were then presented including 
studies of the dissipative anomaly, bottleneck effect and 
thermalisation phenomenon. 

The preliminary results presented here by way of moti- 
vation for the study of the isotropic 3-wave kinetic equa- 
tion have suggested that some further investigation of 
cut-off related phenomena such as the bottleneck would 
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Appendix: Derivation of the 
Aggregation— Fragmentation Equations 

We consider isotropic wave turbulence. The wave 
spectrum , iiy^, is therefore a function of k = jkj only 
which we shall denote by rife. There exist a strightfor- 
ward set of changes of variables which take advantage 
of this isotropy to convert the collision integral over the 
pair of d-dimensional wave- vectors, ki and k2, into one- 
dimensional integrals over frequencies. 

The total number of waves in the system is 

N= rikdk. 

Jm'' 

Denoting the radial coordinate of k G M'' and integration 
over the corresponding angular variables by d&'^^ we can 
write this integral in spherical polar co-ordinates as: 

= n^'^^ / Hkk'^-^dk 

JO 

where we have used the isotropy of nk to integrate 
over the angular variables and denoted the resulting d- 
dimensional solid angle by We now use the isotropy 
of the dispersion relation, 

oj{k)=ck" 

to change variables from integration over k to integration 
over w. 

N = r nuk^-^^dw 

Jo dw 

= c / riojLO " dco. (89) 
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where = 71^(1^). Based on these manipulations, we 
define the angle-averaged frequency spectrum, N^^ by 



= c <i a; a 



(90) 



The angle-averaged frequency spectrum has the advan- 
tage that the total number of waves. A'', and total wave 
energy, E, are given very simply as 



N = [ N^duj 
Jo 

POO 

Jo 



(91) 
(92) 



Our objective is to integrate over angular variables and 
express the kinetic equation entirely in terms of N^. 

We begin by using Eq. (5), to write an evolution equa- 
tion for the time evolution of N^^ : 

f°° dN 

{S4N^]+S2[N^] + S3[N^])cLoi. 



f 

Jo 



This yields a kinetic equation of the form 



dt 



Si[N^]+S2[N^] + S3[N^]. (93) 



To determine the form of Si[N^], 5*2 [A^w] and /SsfA^a;] we 
split the collision integral, ^[nk], into three terms, as 
follows: 



S[nk]dki = aiTi + 02T2 + asTa 



(94) 



(5(a;ki - - Wka) ^(ki - k2 - ka) dki(ik2fik3 

- 47r y"ik2k3ki(«i'^kink3 - a2nk2nk3 - asnw^nw^) 

6{u)k^ - Wks - Wki) <^(k2 - ks - ki) dkidk2dk3 

- 4-77 y ik3kik2(«i'^kink2 - a2nkink3 + a3'^k2»^k3) 

^(a;k3 - Wki - a;k2) <^(k3 - ki - k2) (ikicflc2cflc3 

The variables oi = 02 = ^3 = 1 have been introduced 
simply to indicate which terms are to be grouped to- 
gether. 

Let us first consider the terms proportional to ai. We 
again introduce spherical polar coordinates: 

J dkidk2dks = J {kik2hf-'de\^Ue'^^de'^^dki23, 

(95) 

where the notation dki2^ represents the integration mea- 
sure, dkidk2dk3, over the radial variables. Noting that 



the angular variables only enter into the interaction coef- 
ficients, Lkik2k3 and the k delta functions, we can define 
an angle-averaged interaction coefficient. 



L{k„k2,k3) = 4nJ L^^k2k3'^(ki-k2-k3) (96) 
dQ'i^'>dQ\^^dQ^^\ 

ki k2 K3 ' 

which is a function of the radial variables only. We can 
then write 

T^l = y" £felfe2/C3"fe2"/C3 (^1^2^:3)'^"^ 5(^23) *123 

- j Lk^kski nkinks{kik2k3)^~'^ 5{ujIi) dki23 

- / £fe3feife2 "-fci»^fc2(^lfc2fc3)'^"^ <5(wj2) rffcl23- 



Here the notation S{LOji^) is a compact representation of 
the frequency delta function: 

S{'^}k)=H^i-^j-^k)- (97) 

Now replace the integration over fc's with integration over 

frequencies as we did in Eq. (89) and use Eq. (90) to 
express in terms of N^. The result is 

- Cd y Z^2^3^ilJ2" N^-^N^^5{ujl-^)duji23 (98) 

- Cd y L^^^^^^w^ N^^^N^^J{ujl2)dwi23, 



where 



- ^(^)i(^)i(^)< 



Cd 



(99) 



Finally define 

Ki {wi,ujj) = CdL^i+ojj ojiojj (wi + i^j) ^ • (100) 

Now use the frequency delta-functions to write Eq (98) 
as 

Ti = j Ki{i02,UJ3)N^2N^Ji0J23)duJl23 

Ki{u3,LUi) N^.N^Jiul^) doJi23 (101) 

Ki{ui,U2) N^iN^J{Ljf2) duJi23, 

Comparing with Eq. (93) we see that we should write 

Si [N^]^[Ki (W2, W3) N^^N^J{u!i-LU2-i^3) duJ23 

-lKi{u3,ui) N^,N^,S{u)2-iJ3-coi) rfa;2^102) 

{lVi,U)2) N^iNu2S{ui3-U)i-CJ2) dhJ23- 
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The same set of manipulations can now be applied to the 
terms proportional to a2 and in Eq. (94) to deduce 
the appropriate forms of S2[N^] and 53 [-/V^^]. The results 
are as follows: 



1^23 



S2[Noj]=jK2 {U)2 , W3 ) -/Va^i S{wi-U)2-U}3) dW2 



and 



Ss[Nj\=jK:i {uj2 ,u)3) N^^N^J{u)i-u}2-u}3) du)2 

-jK3{u)3,u)i) N^^N^J{uj2-iO3-0Ji) (ia;2^104) 

,UJ2) N^^N^^5{uJ3-iJi-UJ2) dW23- 

Note that there is a small price to be paid for hiding all 

dependence on d and a: the interaction cocfRcicnts for 
the three collision integrals are, in general, not identical: 

a-d 

' LUi + UJj ' 



. (105) 



Note also that K2{uii,uij) and K3{uJi,u}j), are not sym- 
metric in their arguments although this latter deficiency 
can be removed if desired by symmetrisation. These 
problems are immaterial since the original variables can 
be easily restored if needs be. It is worth noting that 
in the case where d = a, the distinctions between the 
interaction cociEc-ients disappear. Furthermore, even in 
the general case, all three interaction coefficients have the 
same degree of homogeneity. Prom Eq. (100), Eq. (99) 
and Eq. (96), it is easy to work backwards and establish 
that this degree of homogeneity, which we denote by A, 
is 

2/3 -a 



A: 



(106) 



Thus, if one is interested in scaling properties of 3-wave 
kinetic equation, almost everything is determined by this 
single scaling parameter, A. 



Appendix: Transfer integrals 

In this appendix we state the explicit expressions for 
the various transfer integrals which go into the estima- 
tion of the collision integral according to Eq. (80) and 
Eq. (81). 

• Case B : 



Ki{flj,LOk) nk{oJk) dwk 



"of = CljNj I K2{Clj,u)k)ni{Clj+Uk) 



(3) f^" 

^Gy = iljAuj K3{flj,u;k)nk{i0k)ni{flj +u)k) 
■^G^^^ = Aujj Uk K3{flj,ujk)nk{oJk)ni{flj + LOk) 



• Case C : 



AT 

'iL^^ = QjNj / Ki{Q.j,LOk)nk{ijJk)dwk 
= Nj / LOkKi{Qj,LJk)nk{LOk)dwk 
K2{nj,uJk)ni^{nj+oJk) 
£G^^^ = Nj / u)kK2{il.j,u)k)ni^{i}j +LUk) 



lof^ = i}j Aojj I K3{Q.j,ujk)nk{<jJk)ni^{Q.j + ujk) 
fcf^ = Awj I (JkK3{nj,(Jk)nk{uJk)ni^{nj +ujk) 



Cr(l) 

R^j 


= fljNj / Ki{flj,uJk)nk{ii)k)dwk 
Jn* 




Ct(1) 
R^k 


f^k 

= Nj / u)kKi{flj,u!k)nk{LOk) dwk 




R^j 


= flj Nj / K2 {flj , Wfe) riij^ {flj + ujk) 
Jn* 






r<^k 

= Nj / u)kK2{Clj,u)k)ni^{Clj +uik) 




CMS) 
R^o 


= flj Aujj / K3 {^j ,ujk)nk {u!k )ni^{Q.j- 


\-oJk) 


CMS) 
R'-'k 


f^k 

= AiOj / uJkK3{flj,uJk)nk{uJk)nij^{flj - 
Jn* 


\-OJk) 



